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INTRODUCTION 


In  order  to  reduce  the  emphasis  which  is  placed  on  flame 
arresting  tests,  another  reliable  method  for  characterizing 
explosion  suppressants  was  needed.  Therefore,  this  study  was 
undertaken  to  establish  if  a  correlation  could  be  developed 
between  the  air  pressure  drop  characteristics  and  flame  arrestor 
characteristics  for  explosion  suppression  materials.  A  reliable 
test  procedure  was  developed  for  measu- ing  air  pressure  drop 
across  a  foam  sample  in  the  pressure  drop  rig  and  is  presented  in 
Reference  1.  The  pressure  drop  rig  is  shown  in  Figure  1. 

The  flame  tube  rig  is  shown  in  Figure  2.  The  flame  tube 
test  chamber  consists  of  a  rectangular  stainless  steel  tank 
capable  of  containing  combustion  overpressures  as  high  as  120 
psig.  The  test  section  is  square  (1  ft  by  1  ft)  and  is  7.5  ft 
long. 

In  qualifying  a  foam,  tests  are  run  in  accordance  with 
Military  Specification  MIL-B-SSOS^B  (USAF)  in  the  flame  tube 
using  a  propane-air  mixture  at  initial  pressures  of  14.7  and 
17.7  psia.  Examples  of  these  qualifying  runs  are  given  in 
Reference  2. 

Two  types  of  tests  are  run  in  the  flame  tube:  flame 
arrestor  tests  and  explosion  suppression  tests.  Currently  a  fine 
pore  (small  hole)  foam  is  identified  as  a  flame  arrestor  and  the 
coarse  pore  (large  hole)  foam  is  a  combustion  overpressure 
suppressor.  Both  types  of  foam  will  suppress  a  combustion 
overpressure  but  the  fine  pore  foam  will  arrest  the  flame  while 
the  coarse  pore  foam  may  let  the  flame  pass  through. 

In  order  to  develop  a  correlation  between  the  pressure 
drop  tests,  the  results  of  which  can  be  related  to  foam  pore 
size,  and  the  flame  tube  tests,  the  results  of  which  are  also 
related  to  foam  pore  size  through  increased  heat  transfer,  a 
theory  for  the  pressure  rise  in  the  flame  tube  was  developed 
which  includes  heat  loss. 


sure  drop  rig  located  at  Wright-Patterson  Air  Force  Base. 
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Since  it  has  been  established  that  the  performance  of 
arrestors  are  essentially  independent  of  the  material  properties 
of  the  arrestor  and  only  depend  on  cell  dimensions  (see  Reference 
3) ,  it  was  concluded  that  the  theory  could  be  developed 
independent  of  the  material  from  which  the  flame  arrestor  is 
made.  The  theory  does  depend  on  the  fuel-air  mixture  and  the 
initial  temperature  and  pressure.  In  the  next  section  the 
pressure  rise  theory  is  developed  for  the  flame  tube. 

PRESSURE  RISE  THEORY 

The  initial  form  of  the  theory  was  developed  to  closely 
match  the  pressure  rise  that  occurs  in  the  flame  tube  when  a 
mixture  of  propane  and  air  is  burnt  without  a  flame  arrestor 
material  present.  Figure  3  shows  a  schematic  of  the  flame  tube 
showing  the  two  regions  on  either  side  of  the  flame:  the 
products  of  combustion  on  the  left  and  the  unburnt  reactants  on 
the  right.  The  following  simplifying  assumptions  are  made: 

1.  The  gas  is  uniformly  mixed  in  each  region. 

2.  The  pressure  is  constant  throughout  the  tube  at  any 
instant  (i.e.,  the  speed  of  sound  is  much  faster  than 
the  flame  speed,  Vf) . 

3.  The  reactants  are  compressed  in  a  reversible  process. 

4.  The  number  of  moles  of  products  is  approximately  equal 
to  the  number  of  moles  of  reactants. 

5.  Ideal  gas  equations  are  valid  for  both  regions. 

6.  The  specific  heats  of  the  reactants  and  products  are 
approximately  equal. 

7.  The  initial  temperature  of  the  gas  is  equal  to  the  wall 
temperature . 

8.  The  wall  temperature  remains  constant. 
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The  conservation  of  mass  for  the  two  regions  of  Figure  3  is 


mp  +  mu  -  m  ( 1  ) 

Since  the  volume  is  also  constant,  we  have  for  the  two  regions 

Vp  ♦  Vu  -  V  (2) 

Using  the  specific  volumes  in  Equation  (2)  and  combining  it 
with  Equation  (1)  yields  after  rearrangement 

v^'1)  ■"“ft”)  (3> 

Since  the  pressure  in  both  regions  is  the  same  and  the  number  of 
moles  are  assumed  to  be  equal  then  it  follows  from  the  ideal  gas 
equation  that 


In  view  of  Equations  (3)  and  (4)  it  follows  that 


(5) 


The  energy  equation,  including  the  heat  released  by  the  reaction, 
is 


mp^ v ( Tp 


T  )  -  mDH  + 
w  P 


m  C  (T 
u  V  u 


-  V  - 


(6) 


We  have  assumed  that  the 
and  the  products  are  the 
related  to  the  adiabatic 


specific  heats  of  the  unburnt  reactants 
same.  The  value  of  H  in  Equation  (6)  is 
flame  temperature,  T^,  by  the  following 


H  =  C  (T.  -  T  ) 
v  f  w 


(7) 
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Substituting  Equations  (7)  and  (1)  into  (6)  and  rearranging 
yields 


VTp  \  ®p 

ra  Wu  /  "  ■ 


<Tf  -  Tw) 


w 


CvmTw 


(8) 


In  view  of  Equation  (5)  above 


m 


C  mT 
v  w 


but  from  the  ideal  gas  equation 


(9) 


(10) 


Hence,  it  follows  from  Equations  (9)  and  (10)  that 


(11  ) 


The  rate  of  product  production  is  related  to  the  flame 
speed  Vf  as 


dmp 

dt~ 


VuAt 


(12) 


According  to  Reference  4,  the  laminar  flame  speed  is  independent 
of  pressure  for  second-order  reactions  which  occur  for  many 
hydrocarbon  fuels.  However,  the  flame  speed  is  a  function  of  the 
temperature,  T  ,  in  the  unburned  reactants.  According  to 
Reference  4,  the  flame  speed  is  given  by 


V 


f 


(13) 


where  the  exponent  mg  ranges  between  1.5  and  2.  Although  the 
turbulent  flame  speed  is  much  greater  than  the  laminar  flame 
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In  view  of  Equation  (11)  for  zero  heat  transfer  we  have 


dt  T  dt 


w 


(19) 


If  we  substitute  Equation  (18)  into  (19)  and  integrate  the 
resultant  equation  the  solution  is 


[■ 


(1  -  me)(Y  -  1) 


me'TTf~Ty 


(20) 


Equation  (11),  for  zero  heat  transfer,  can  be  used  to  solve  for 
the  mass  of  products 


m 


P/P„  -  1 
o 

T  7t  - 
1f/  xw 


(21  ) 


Equation  (21 )  shows  that  if  mD/m  -  1 ,  then  P/P^  -  Tp/T 

r  O  I  W 

which  it  should  in  an  adiabatic  constant  volume  burning  of  the 
fuel.  Since  the  unburnt  reactants  are  compressed  isentropically 
it  follows  that 


In  view  of  Equation  (5)  we  have 


+ 


v  /  v 


u 

/m 


1 


(22) 


(23) 


Equations  (20)  to  (23)  along  with  Equation  (17)  give  a 
complete  solution  of  the  rise  time  histories  of  all  the 
variables.  The  equations  were  used  to  fit  the  data  shown  in 
Figure  4. 


r, 
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Figure  4.  Overpressure  versus  time  in  the  flame  tube  without  foam. 


The  curve  marked  "1"  in  Figure  4  was  plotted  from  data 
obtained  in  a  test  conducted  during  the  current  program.  The 
curve  marked  w2"  in  Figure  4  was  plotted  from  data  presented  in 
Reference  5.  The  higher  peak  of  curve  "1"  of  Figure  4  is  a 
result  both  of  the  higher  initial  pressure  used  in  that  test  and 
the  higher  volume  percentage  of  propane  in  air.  (As  shown  in 
Reference  5,  the  maximum  overpressure  occurs  at  5  percent  by 
volume  for  propane  in  air.) 

In  order  to  fit  the  data  shown  in  Figure  4  the  following 
procedure  was  used. 

1.  A  value  of  the  exponent,  mg,  of  Equation  (13)  was 
chosen  between  1.5  and  2. 

2.  The  time  of  the  occurrence  of  the  peak,  tp ,  is  obtained 
from  the  data  (e.g.,  tp  -  0.44  for  curve  2  of  Figure 
4). 

3.  The  ratio  of  the  peak  pressure  to  the  initial  pressure 
(P/PQ)  in  an  adiabatic  constant  volume  burning  is  equal 
to  the  ratio  of  the  adiabatic  flame  temperature  to  the 
initial  temperature  (T„/T  ). 

I  W 

4.  Using  Equation  (20)  with  procedures  (2)  and  (3)  above 
solves  for  the  initial  value  of  the  flame  speed: 


Using  the  above  procedure  the  data  of  Table  1  were 
calculated.  The  data  shown  in  Table  1  were  calculated  at  0.01  s 
time  intervals  from  0  s  to  0.44  s.  The  last  column  of  Table  1 
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TABLE  1 


TIME  HISTORIES  OF  VARIOUS  PROPERTIES  CALCULATED 
FROM  THE  ANALYTICAL  SOLUTION  WITHOUT  HEAT  LOSSES 


BCPONOir  iff=  1.5  Vfo=  4.341  ft/s  T£=4050.0  deg-R  Tw=  500.0  deg-R 
Po=  17.14  psia  L=  7.5  ft  t-peak=0.44  s 


t 

dP 

P 

Tr 

T*’ 

Rho/RhoO 

0.00 

0.00 

17.14 

500.00 

500.00 

1.000 

0.00 

0.01 

0.72 

17.86 

505.92 

3056.52 

1.030 

0.26 

0.02 

1.48 

18.62 

511.95 

3077.44 

1.061 

0.52 

0.03 

2.27 

19.41 

518.09 

3098.45 

1.093 

0.77 

0.04 

3.10 

20.24 

524.34 

3119.57 

1.126 

1.01 

0.05 

3.97 

21.11 

530.70 

3140.79 

1.161 

1.25 

0.06 

4.89 

22.03 

537.18 

3162.10 

1.196 

1.48 

0.07 

5.85 

22.99 

543.78 

3183.52 

1.233 

1.71 

0.08 

6.86 

24.00 

550.50 

3205.04 

1.272 

1.94 

0.09 

7.92 

25.06 

557.34 

3226.65 

1.312 

2.16 

0.10 

9.04 

26.18 

564.32 

3248.36 

1.353 

2.37 

0.11 

10.21 

27.35 

571.42 

3270.18 

1.396 

2.58 

0.12 

11.44 

28.58 

578.66 

3292.09 

1.441 

2.78 

0.13 

12.74 

29.88 

586.04 

3314.10 

1.487 

2.99 

0.14 

14.10 

31.24 

593.56 

3336.20 

1.535 

3.18 

0.15 

15.54 

32.68 

601.23 

3358.41 

1.586 

3.37 

0.16 
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gives  the  position,  Lp ,  of  the  flame  at  the  value  of  time  shown 
in  the  first  column.  The  position  is  calculated  from  the 

equation 


It  is  of  interest  to  plot  the  overpressure  given  in  column 
2  versus  the  relative  flame  position  Lp/L  as  shown  in  Figure  5. 

We  have  referred  to  the  relative  flame  position  as  the  void 
fraction  since  in  the  flame  suppression  tests  the  flame  tube  is 
packed  with  foam  from  some  position  in  the  flame  tube  to  the  end 
of  the  tube.  The  volume  in  front  of  the  foam  over  the  total 
volume  of  the  tube  is  referred  to  as  the  void  fraction.  If  a 
flame  traveled  up  to  the  position  of  the  foam  and  was 
extinguished  at  that  point,  the  value  of  the  overpressure 
achieved  would  be  given  by  the  curve  shown  in  Figure  5.  This 
concept  will  be  compared  to  experimental  results  in  a  later 
section . 

The  overpressure  shown  in  Table  1  is  plotted  versus  time  in 
Figure  6  where  it  is  compared  to  the  data  shown  in  Figure  4.  The 
comparison  is  reasonably  good  for  both  curves.  The  peak  value  is 
higher  than  the  experimental  value  by  about  11  percent  which  is 
probably  due  to  neglecting  the  heat  transfer.  The  effects  of 
heat  transfer  will  be  discussed  in  a  later  section.  However,  in 
the  next  section  we  will  first  develop  a  theory  for  the  flame 
arrestor  tests  and  compare  the  results  to  experiments. 

FLAME  ARRESTOR  TESTS 

Figure  7  shows  the  configuration  of  the  flame  tube  for 
flame  arrestor  testing.  The  front  of  the  foam  is  located  15 
inches  from  the  front  of  the  flame  tube  which  gives  a  combustion 
void  of  16.7  percent.  The  region  behind  the  foam  is  the  relief 
volume.  The  thickness  of  the  foam  used  in  the  test  varies  from  5 
inches  down  to  0.5  inch  or  the  smallest  thickness  of  foam  that 


r 


does  not  blow  through.  During  the  burning  the  flame  front  drives 
the  combustion  void  into  two  regions:  one  contains  the  products 
of  combustion  and  the  other  contains  the  unburnt  mixture. 

Two  extreme  effects  of  the  foam  can  be  considered.  The 
first  extreme  case  is  represented  by  a  foam  that  does  not  prevent 
the  transmission  of  the  pressure  waves  through  the  foam  but 
extinguishes  the  flame  when  it  reaches  the  foam.  In  this  case 
the  overpressure  would  be  about  4  psid  which  can  be  obtained  from 
Figure  5  at  a  void  fraction  of  1 6 . 7  percent. 

In  the  other  extreme  case  all  of  the  pressure  waves  would 
be  reflected  by  the  foam  (it  would  act  like  a  solid  wall)  and  the 
pressure  would  reach  the  peak  value  consistent  with  adiabatic 
constant  volume  burning. 

It  was  believed  that  the  foam  acted  somewhere  between  these 
two  extremes  and  would  show  a  resistance  that  would  increase  the 
overpressure  in  front  of  the  foam  above  the  value  shown  at  16.7 
percent  void  fraction  in  Figure  5.  A  theory  was  developed  for 
this  case,  which  is  similar  to  one  developed  in  Reference  6,  as 
follows.  However  the  experimental  results,  as  discussed  later, 
show  that  foam  does  not  substantially  attenuate  pressure  waves 
and  therefore,  the  first  extreme  case  is  essentially  what 
occurs . 

The  conservation  of  mass  for  a  control  volume  in  front  of 
the  foam  leads  to  an  equation  similar  to  Equation  (5): 


(26) 


where  mf  is  the  mass  flow  rate  through  the  foam  and  the  integral 
gives  the  total  mass  flow  out  of  the  region  in  front  of  the  foam 
during  the  time  t. 


When  the  energy  equation  for  the  control  volume  is  combined 
with  Equation  (26),  the  resulting  expression  for  the  pressure  is 


It  can  be  shown  that  the  flow  from  a  reservoir  through  a  pipe 
with  and  without  foam  in  the  pipe  leads  to  tne  equation 


where  m^  is  the  mass  flow  rate  for  an  inviscid  case  (i.e.,  no 
foam)  and  fT0T  is  the  total  friction  factor  of  the  foam.  The 
mass  flow  rate  at  any  location  in  the  flame  tube  without  foam  can 
be  related  to  the  time  rate  of  change  of  the  pressure  as  follows. 
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(29) 


Combining  Equations  (27),  (28),  and  (29)  leads  to 


(30) 


Equation  (30)  reduces  to  Equation  (19)  if  both  LR  and  QT  are 


The  pressure  behind  the  foam  is  related  to  the  mass  flow 
into  the  region  R  by  Equation  (29)  if  P  is  replaced  by  PR  and  mfi 
is  replaced  by  m^..  It  follows  that 
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TOT 


(31  ) 


Since  Equation  (18)  is  still  valid,  Equations  (30)  and  (31) 
can  be  solved  for  zero  heat  transfer.  Results  from  the  solution 
are  shown  in  Figure  8.  The  theoretical  pressure  rise  in  front  of 
the  foam  would  depend  on  the  thickness  of  the  foam  since  fTQT  in 
the  equations  above  is  directly  proportional  to  the  foam  thick¬ 
ness.  As  seen  from  Figure  8,  the  pressure  in  front  of  the  foam 
is  predicted  to  be  substantially  higher  than  the  pressure  behind 
the  foam  when  we  use  values  of  f  found  with  the  pressure  drop 
rig. 

We  first  compared  the  theory  to  data  obtained  from 
Reference  2.  The  data  points  plotted  between  0.5  and  2  inches 
in  Figure  8  are  from  Reference  2  for  Crest  type  III  red  foam. 
Although  two  pressure  transducers  were  used  in  the  tests,  they 
were  located  behind  the  foam  and  therefore,  the  results  were  not 
contradictory  to  the  theory. 

In  order  to  conduct  definitive  tests  we  ran  tests  with 
3  and  5  inch  thick  pieces  of  foam.  The  pressure  transducers  were 
relocated  so  that  one  was  in  front  of  the  foam  and  the  other  was 
behind  the  foam.  The  results  of  these  tests  are  also  plotted  in 
Figure  8. 

These  tests  demonstrated  convincingly  that  the  foam  does 
not  hinder  the  passage  of  pressure  waves  through  it  since  both 
transducers  had  essentially  the  same  time  histories.  Therefore 
the  peak  pressure,  provided  the  flame  is  extinguished  by  the 
foam,  can  be  predicted  reasonably  well  from  the  location  of  the 
foam  using  Figure  5. 
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It  should  be  pointed  out  that  the  minimum  foam  thickness 
that  will  extinguish  the  flame  is  fairly  thin.  Therefore,  the 
method  used  to  support  the  foam  in  these  tests  has  a  significant 
effect  on  the  results.  Tests  reported  in  Reference  2  showed 
erratic  results  until  it  was  determined  that  one  of  two  racks 
being  used  gave  inadequate  support  of  the  foam.  After  this  rack 
was  discarded,  the  results  were  consistent  and  foam  extinguished 
the  flame  even  when  tested  down  to  a  one-half  inch  thickness. 

In  the  next  section  results  of  the  explosion  suppressant 
testing  will  be  compared  to  the  theory. 

EXPLOSION  SUPPRESSANT  TESTS 

Figure  9  shows  a  schematic  of  the  flame  tube  configuration 
used  for  explosion  suppressant  tests.  In  these  tests  various 
combustion  voids  are  tested  by  packing  various  amounts  of  foam 
into  the  flame  tube. 

Results  from  a  number  of  tests  described  in  Reference  2  are 
shown  plotted  in  Figure  10  where  the  peak  overpressure  is  plotted 
versus  void  fraction.  The  data  are  compared  to  the  theoretical 
curve  shown  in  Figure  5. 

In  Figure  10  two  test  series  of  Crest  type  II  red  foam  are 
shown.  In  the  first  series  tests  were  run  with  racks  present  in 
the  void  space  in  front  of  the  foam.  As  seen  in  Figure  10,  the 
results  were  very  erratic  with  some  very  high  overpressures.  The 
data  for  the  other  two  types  of  foam  were  also  obtained  with  the 
racks  present  but  each  of  these  foams  was  qualified  under 
Military  Specification  MIL-B-83054B  (USAF).  Therefore  these  two 
types  were  not  rerun  without  racks. 

After  the  initial  tests  for  the  type  III  foam  were  run  with 
the  racks,  it  was  determined  that  earlier  qualification  tests 
were  run  on  other  foams  without  the  racks  present.  Apparently 
the  racks  moved  during  the  test  and  compressed  the  foam  which 
enlarged  the  combustion  void  and  consequently  increased  the 


Configuration  of  flame  tube  for  explosion 
suppressant  testing. 


overpressure.  Subsequently  the  qualification  tests  for  the  type 
III  foam  were  rerun  without  the  racks  in  the  combustion  void. 

As  can  be  seen  in  Figure  10,  the  results  were  no  longer 
erratic  and  they  agreed  very  well  with  the  theory  developed 
earlier.  Because  of  these  results  it  was  decided  to  rerun  one  of 
the  other  foams  without  racks  to  see  if  the  results  would  more 
nearly  agree  with  the  theory.  The  type  I  orange  foam  was  chosen 
for  testing  without  racks  and  the  results  are  shown  in  Figure  11 
where  both  the  type  III  and  type  I  results  without  racks  are 
shown . 

At  combustion  voids  of  50  percent  or  less  the  results  are 
in  excellent  agreement  with  the  theory.  Above  50  percent  the  one 
point  for  the  type  III  red  foam  still  agrees  quite  well  with  the 
theoretical  curve  while  the  two  points  shown  for  the  type  I  foam 
show  a  higher  peak  pressure  than  the  theory  would  predict.  This 
may  be  because  the  foam  did  not  extinguish  the  flame  and  burning 
continued  into  the  foam.  This  phenomena  is  not  uncommon  in  the 
explosion  suppressant  testing  and  was  observed  on  a  number  of 
occasions  during  testing. 

The  type  I  orange  foam  is  the  most  open  foam  with  a  nominal 
10  pores  per  inch  and  was  chosen  for  retesting  because  of  this. 
The  type  III  red  foam  has  a  nominal  25  pores  per  inch  and 
therefore,  has  a  smaller  pore  size  than  the  type  I  orange  foam. 
The  type  I  orange  foam  would  allow  the  flame  to  penetrate  the 
foam  more  readily  than  the  type  III  red  foam. 

The  larger  combustion  voids  bring  the  mixture  within  the 
foam  to  a  higher  pressure  and  temperature  before  the  flame 
reaches  the  foam.  This  preheating  may  permit  the  flame  to 
penetrate  into  the  foam  and  increase  the  peak  pressures  in  the 
tests  at  the  higher  values  of  combustion  void.  The  data  for  the 
type  I  orange  foam  are  commensurate  with  this  concept.  The  type 
III  red  foam,  because  of  its  smaller  pores,  would  require  more 
preheating  than  the  type  I  orange  foam.  This  may  explain  why  the 
type  III  red  foam  still  agrees  with  the  theory  at  a  combustion 


Figure  11.  Comparison  of  experimental  data  tc  theory  for  two  types  of  foam 
Test  data  taken  without  racks  present  in  the  void  space. 
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void  of  60  percent  while  the  type  I  orange  foam  exhibits  an 
increase  in  overpressure  over  the  theory  at  a  combustion  void  of 
60  percent  as  shown  in  Figure  11. 


In  order  to  extend  the  theory  to  investigate  the  propaga¬ 
tion  of  the  flame  into  the  foam  it  is  necessary  to  consider  the 
effects  of  heat  transfer. 


HEAT  LOSS  STUDY 


When  the  effects  of  heat  transfer  are  considered,  the  heat 
loss  terms  in  Equations  (11)  and  (16)  must  be  retained.  A  closed 
form  solution  was  not  found  for  this  case,  but  a  fourth-order 
Runge-Kutta  numerical  procedure  was  used  to  solve  Equations  (14), 
(16),  and  (19)  (a  heat  transfer  term  was  added  to  Equation  (19)). 
The  method  was  checked  with  the  exact  solution  by  setting  the 
heat  transfer  to  zero.  The  heat  transfer  term,  required  in 
Equation  (19),  is  obtained  from  Equation  (11)  by  differentiating 
with  respect  to  time.  The  result  is 


<0  Tf  -  T„  d0r ) 


dt 


w 


dt 


The  heat  transfer  term  in  Equation  (32)  is  the  total  value 
and  includes  the  heat  loss  from  the  products  and  from  the  unburnt 


reactants  ($u  of  Equation  (16)).  The  standard  thermodynamics 


convention  was  used  for  the  heat  transfer  terms:  heat  into  the 
control  volume  is  positive.  The  heat  transfer  from  the  unburnt 
reactants  is 


C  mT 
u  w 
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Hence,  the  term  needed  for  Equation  (16)  is 
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The  density  ratio  can  be  eliminated  in  favor  of  the 
temperature  and  pressure  ratios  resulting  in 
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(35) 


In  a  similar  way  the  heat  transfer  term  needed  for  Equation  (32) 
can  be  obtained 
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The  terms  in  square  brackets  contain  the  heat  transfer  coeffi¬ 
cients  for  each  region  and  a  geometric  parameter  (C.D/A.)  which 

In  u 


is  related  to  the  hydraulic  diameter  (D 
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Mt/CiR). 


Thus,  for 


example 
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which  shows  the  significance  of  the  hydraulic  diameter  and 
consequently  the  small  pore  sizes  in  increasing  the  heat 


transfer.  For  small  enough  values  of  Dh  the  effect  of  the 


cooling  is  so  large  that  the  flame  is  quenched.  If  the  flame  is 


not  quenched,  the  effect  of  the  cooling  for  small  is  to  reduce 


the  overpressure  substantially  as  occurs  in  the  explosion  sup¬ 
pression  foams  when  flame  propagates  through  the  foam. 

* 


In  the  above  equations  the  heat  transfer  coefficient,  h  , 
is  a  function  of  temperature  and  can  be  a  result  of  both  free 


convection  and  radiation.  For  gray  body  radiation  the  value  of 
h*  in  BTU/hr/f t2/°R  is 
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The  results  of  using  Equation  (38)  with  various  values  of 
the  emissivity,  e,  are  shown  in  Figures  12  through  16  where  the 
theory  is  compared  to  the  experimental  data.  The  page  parameter 
is  the  initial  flame  speed  which  varies  from  4  ft/s  for  Figure  12 
to  6  ft/s  for  Figure  16  in  steps  of  0.5  ft/s.  The  best  fits  are 
obtained  with  a  flame  speed  of  about  5  ft/s  and  e  of  about  0.3. 
However,  when  compared  to  the  data,  the  rate  of  cooling  appears 
to  be  too  high  during  the  pressure  rise  phase  and  not  high  enough 
during  the  cooling  phase. 

In  an  attempt  to  obtain  a  better  fit  to  the  data  the  gas 
radiation  theory  developed  by  Hottel  and  Egbert  (References  7.  8» 
and  9)  was  used.  They  present  a  series  of  charts  to  account  for 
the  radiation  due  to  water  vapor  and  carbon  dioxide.  These 
charts  are  widely  used  to  determine  the  emissivity  of  gases. 

The  following  equations,  which  approximate  the  charts,  were 
given  by  Ganapathy  in  Reference  10  to  approximate  e. 

-KL. 

e  -  0.9(1  -  e  D)  (39) 

where 


Lb  -  3.4  AtL/L-CiR  =  3.4  Dh/4  (41) 

In  Equations  (39)  through  (41)  T  is  the  temperature  in 
Kelvin,  Pc  and  Py  are  the  partial  pressures  of  carbon  dioxide  and 
water  vapor  in  atmospheres,  and  is  the  beam  length  in  meters. 

A  term  was  also  added  to  account  for  free  convection. 
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e  with  gray  body  radiation  curves  for  var 
flame  speed  is  4  ft/s. 
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Figure  15-  Pressure  versus  time  with  gray  body  radiation  curves  for  various  value 


Using  the  well  known  Sutherland  formula  for  p  and  k  and  recalling 
that  8  is  1 /T  for  an  ideal  gas  yields 


(43) 


Using  a  value  of  C  -  0.55  and  substituting  for  the  various  values 

* 

yields  a  value  of  1.3  for  the  term  in  square  brackets  if  h~c  is 
in  BTU/hr/f t2/°R. 

The  results  of  the  above  approach  are  seen  in  Figure  17 

where  the  theory  is  compared  to  experiment.  The  fit  shown  was 

for  a  flame  speed  of  4.34  ft/s.  The  cooling  does  not  follow  the 

* 

data  very  well  and  apparently  the  value  of  h  is  too  small. 

The  gas  radiation  theory  used  was  originally  substantiated 
by  experiments  at  high  gas  temperatures.  Therefore,  it  is 
believed  that  the  heat  transfer  from  convection  must  be 
substantially  higher  than  that  predicted  from  Equation  (43). 

Since  a  purely  theoretical  approach  did  not  yield  satisfactory 
results,  an  emperical  fit  was  obtained  using  the  equation 

h*  -  C ( T  -  T  )m  (44) 

w 

During  the  cool  down  phase  the  temperature  of  the  products 
can  be  determined  from  the  following  equation. 


Substitution  of  Equation  (44)  in  Equation  (45)  yields 
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In  view  of  Equation  (44)  we  can  eliminate  C  in  favor  of  the 

» 

initial  value  of  heat  transfer  coefficient  h^ . 


dT 

dt 


(47) 


The  solution  to  Equation  (47)  is 
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(48) 


By  trial  and  error  a  fit  of  Equation  (48)  was  made  to  the 

cool  down  data  and  is  also  shown  in  Figure  17.  The  fit  to  the 

* 

data  was  excellent  for  a  value  of  m  =  1.2795  and  Ik  *=  33.39 
BTU/hr/f t2/ °R  (or  C  -  0.001129).  The  value  of  P^/Po  as  read  from 
the  data  was  7.24. 

Although  time  was  not  available  under  the  current  contract, 
Equation  (44)  should  be  used  along  with  the  hydraulic  diameter  of 
the  various  foams  to  study  the  effects  of  burning  within  the  foam 
and  for  determining  when  the  flame  will  be  extinguished  by  the 
foam.  The  relationship  for  determining  the  pore  size  required 
for  extinguishing  the  flame  will  be  discussed  in  the  next 
section . 

PORE  SIZE  CORRELATION  WITH  PRESSURE  DROP 

In  determining  if  a  foam  will  extinguish  a  flame  or  if  a 
flame  will  propagate  into  the  foam  it  is  necessary  to  determine 
the  pore  size  of  the  foam. 


A  criteria  for  extinguishing  a  flame  was  presented  by 
Wilson  and  Attalah  in  Reference  11.  Flame  quenching  is  related 


to  the  hydraulic  diameter  as  would  be  anticipated  from  the  heat 

loss  study  presented  above.  This  quenching  only  occurs  if  the 

* 

hydraulic  diameter  is  less  than  a  critical  value,  Dh. 

* 

Wilson  and  Attalah  relate  to  the  thermal  diffusivity  a, 
laminar  flame  speed  ,  and  the  Peclet  number,  Npg: 


*  a 

Dh  ‘  NPe 


(49) 


According  to  Wilson  and  Attalah  for  hydrocarbon-air 
mixtures  at  standard  conditions  a  typical  value  of  Peclet  number 
is  60.  Therefore  the  critical  hydraulic  diameter  for  many  fuels 


Reference  4  develops  a  similar  expression  for  .  The  idea 
of  a  minimum  diameter  below  which  a  fuel-air  mixture  cannot 
support  a  flame  is  so  widely  accepted  that  the  "quenching 
diameter"  is  often  listed  with  the  flame  speed  as  a  fundamental 
property  of  the  fuel.  Both  Reference  11  and  4  give  values  of 
0.25  cm  (0.1  in)  as  typical  values  of  hydrocarbon  fuel.  (Note 
that  Reference  11  suggests  that  distances  cited  in  the  literature 
for  parallel  plates  or  slots  be  multiplied  by  1.4  to  obtain  the 
hydraulic  diameter.) 

* 

For  propane  Reference  11  gives  a  value  of  0.2J|  cm  for  D^. 
They  also  suggest  for  fuel-air  mixtures  having  laminar  flame 
speeds  up  to  100  cm/s  (3.3  ft/s)  that  a  conservative  design 
criterion  is 

D*  <  0.075  cm  (0.03  in) 


which  has  a  safety  factor  of  two  or  more,  depending  on  the  fuel 
(hydrogen  and  acetylene  may  require  smaller  values  of  ) . 


* 

The  value  of  depends  on  the  pressure  and  the  temperature 

of  the  unburnt  reactants.  For  hydrocarbon  fuels  burning  in  air 

* 

the  following  adjustments  to  Dh  should  be  made: 


For  propane  the  value  of  is  0.24  cm  for  PQ  equal  to  1 

o 

atmosphere  and  T,,_  at  a  temperature  of  80°F.  Equation  (51)  can 

uo  # 

be  used  to  adjust  to  other  conditions. 

The  hydraulic  diameter  of  the  foam  can  be  estimated  from 
the  formula 


_ 4- volume 

surface  area 


4  Vol  4 

"*7"  '  s 


(52) 


The  surface  area  to  volume  ratio,  s,  is  the  specific 
surface  area  of  the  foam.  For  foam  with  a  density  of  1.75  pounds 
per  cubic  foot  Scott  (form  3645-R-0375)  shows  the  variation  of 

the  specific  surface  area  as  a  function  of  the  pores  per  inch.  A 

satisfactory  fit  to  their  data  is 

S  -  9.64  Ppi1 *169  (53) 

2  3 

where  S  has  units  of  ft  /ft  and  P  is  in  pores  per  inch.  Thus, 
in  view  of  Equation  (52)  an  estimate  of  the  hydraulic  diameter  is 

D.  -  0. 41  49  P  <_1  ' 1  69  (54) 

h  pi 

where  is  in  feet.  If  the  constant  0.4149  is  replaced  by 

1.960,  the  result  would  be  in  centimeters.  If  Equation  (54)  is 

used  to  calculate  the  hydraulic  diameter  of  the  type  I  orange 

foam  (10  pores  per  inch)  and  the  type  III  red  foam  (25  pores  per 

inch),  the  result  is  0.133  cm  and  0.0455  cm,  respectively.  Thus, 

both  foams  have  hydraulic  diameters  less  than  the  0.24  cm  value 
*  * 
for  Dh  at  1  atmosphere  and  80°F.  However,  if  the  value  of  is 
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adjusted  by  Equation  (51)  using  the  data  of  Table  1,  then  the 

result  can  be  plotted  as  shown  in  Figure  18.  Figure  18  shows  the 

# 

variation  in  in  centimeters  as  the  flame  progresses  down  the 

flame  tube.  When  the  flame  has  traveled  to  4.5  feet  (60  percent 

* 

of  7.5),  the  new  value  of  is  0.075  cm.  Thus,  the  flame  would 

be  expected  to  penetrate  the  type  I  orange  foam  but  not  the  type 

III  red  foam  just  as  the  experiments  indicated  (see  Figure  11). 

« 

At  the  50  percent  location,  from  Figure  18,  is  0.092  cm.  This 
would  indicate  that  the  type  I  orange  foam  would  also  burn 
through  at  50  percent  but  the  data  of  Figure  11  indicate  that  it 
did  not. 

An  empirical  fit  was  made  to  the  pressure  drop  data  by 
using  the  information  in  the  Scott  circular  (form  3645-R-0375) 
with  the  following  result. 

f  -  5.352  Pp.0*8481  5  (55) 

The  values  of  the  constants  were  determined  using  average  values 
of  f  at  20  and  45  pores  per  inch.  Solving  for  Ppi  in  terms  of 
f COR  Pv/P  ”  59021  per  ft)  yields 

Ppl  -  0.965T  fCoRK179  <56) 


Table  2  shows  calculated  values  of  Ppi  from  Equation  (56) 
for  a  number  of  foam  samples  that  were  cut  from  different 
locations  in  the  bun.  The  nominal  values  of  the  pores  per  inch 
for  each  sample  are  also  shown  in  Table  2.  The  agreement  is 
reasonable  from  15  to  45  pores  per  inch.  Thus,  pressure  drop 
data  can  be  used  to  determine  the  value  of  fCQR  which  can  be  used 
to  determine  the  hydraulic  diameter.  If  Equation  (56)  and  (54) 
are  combined,  the  result  is 


0.4322  f 


-1.372 


COR 


(57) 
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Figure  18.  Critical  value  of  the  hydraulic  diameter  for  propane  versus 
flame  position  in  the  flame  tube. 


TABLE  2 


CALCULATED  PORES  PER  INCH  COMPARED  TO  THE  NOMINAL 
VALUE  FOR  VARIOUS  FOAM  SAMPLES 


FOAM  ID 

NOMINAL  PORES 
PER  INCH 

fCOR 

CALCULATED  PORES 
PER  INCH 

T15-FV 

15 

10.20 

14.9 

-REV 

i 

10.01 

14.6 

-RV 

12.17 

18.4 

-LV 

9.36 

13.5 

-BH 

I 

8.15 

11  .5 

-MH 

.07 

13.0 

-TH 

15 

9.47 

13.7 

T4-A1M1 

15 

8.08 

1  1  .5 

-A1M2 

• 

8.50 

12.2 

-B2M1 

8.26 

11.8 

-B2M2 

8.45 

12.1 

-C3M1 

! 

8.23 

1  1  .7 

-C3M2 

15 

8.62 

12.4 

T20-RV 

20 

13.63 

21  .0 

-LV 

15.93 

25.3 

-BH 

1 

9.25 

13.3 

-DB 

14.08 

22.0 

-MH 

! 

11.54 

17.0 

-TH 

11  .73 

18.0 

-REV 

20 

15.32 

24.0 

T3-A1M1 

25 

15.02 

24.0 

-A1M2 

1 

14.84 

23.7 

-C2M1 

14.69 

23.4 

-C2M2 

14.63 

23.3 

-E3M1 

1 

14.08 

22.2 

-E3M2 

25 

13.85 

21  .8 

T5-A1M1 

25 

16.17 

26.2 

-A1M2 

1 

15.30 

24.5 

-B2M1 

14.37 

22.8 

-B2M2 

15.28 

24.5 

-C3M1 

f 

16.03 

25.9 

-C3M2 

25 

16.08 

26.0 

W742-V1 -9 

45 

24.73 

42.4 

The  value  of  D.  calculated  from  Equation  (57)  can  be  compared  to 

*  “ 

determined  with  Figure  18  to  estimate  at  what  void  fraction 
the  foam  would  no  longer  extinguish  the  flame.  For  the  cases 
when  the  foam  does  not  extinguish  the  flame,  the  value  of  D^, 
calculated  from  Equation  (57),  can  be  used  with  the  equations 
from  the  heat  loss  studies  to  determine  the  effect  on  the  final 
value  of  the  overpressure. 

RESULTS  AND  RECOMMENDATIONS 

A  number  of  useful  results  were  obtained  during  this 
research  effort: 

•  A  correlation  was  made  between  the  number  of  pores  per 
inch  and  hydraulic  diameter  of  the  foam. 

•  A  correlation  was  made  between  the  friction  coefficient 
f Cop  and  the  number  of  pores  per  inch  of  the  foam. 

•  A  generalized  flame  quenching  criteria  was  presented. 

•  The  variation  of  the  critical  hydraulic  diameter  with 
pressure  and  temperature  was  presented. 

•  A  procedure  for  relating  the  critical  hydraulic  diameter 
to  the  flame  position  in  the  flame  tube  wa3  developed. 

•  The  void  fraction  required  to  produce  burn  through  the 
foam  was  determined. 

•  An  empirical  equation  was  developed  for  the  cooling 
phase  in  the  flame  tube. 

•  A  theory,  that  compared  favorably  to  experiments,  for 
predicting  flame  tube  overpressure  was  developed. 

•  It  was  shown  that  the  flame  arrestor  results  are 
independent  of  fCQR  (or  pore  size)  and  thickness 
provided  it  quenches  the  flame. 

•  Pressure  waves  are  transmitted  through  the  foam  with 
little  or  no  attenuation. 


•  A  pressure  rise  theory  was  developed. 

•  The  equations  and  procedures  presented  enable  the 
performance  of  the  flame  tube  to  be  determined  from  the 
values  of  the  initial  temperature,  initial  pressure, 
void  fraction,  and  hydraulic  diameter  of  the  foam. 

It  is  recommended  that  the  known  variation  of  the  flame 
speed  and  adiabatic  flame  temperature  be  integrated  into  the 
theory.  Also,  the  theory  should  be  extended  to  determine  the 
overpressure  that  results  when  the  flame  continues  to  burn 
through  the  foam. 
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